Signatures of a staggered-fiux phase in the t-J model with two holes on a 32-site lattice 
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We study the relevance of the staggered-fiux phase in the t-J model using a system with two 
holes on a 32-site lattice with periodic boundary conditions. We find a staggered-fiux pattern in 
the current-current correlation in the lowest energy d-wave state where there is mutual attraction 
between the holes. This staggered correlation decays faster with distance when the hole binding 
becomes stronger. This is in complete agreement with a recent study by Ivanov, Lee and Wen 
(Phys. Rev. Lett. 84, 3958, (2000)) based on the 5(7(2) theory, and strongly suggests that the 
staggered-fiux phase is a key ingredient in the t-J model. We further show that this staggered-fiux 
pattern does not exist in a state where the holes repel each other. Correlations of the chirality 
operator Si • (S2 x S3) show that the staggered pattern of the chirality is closely tied to the holes. 



PACS: 7f.27.+a, 71.10.Fd, 75.40.Mg 

The t-J model is one of the most important micro- 
scopic models in the study of high temperature supercon- 
ductivity. It describes a doped Cu02 plane as a system 
of holes (or Zhang- Rice singlets 0]) moving in a spin 
background described by the antiferromagnetic Heisen- 
berg model. The Hamiltonian is 
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where cJ CT and Ci G are the projected fermion operators, 



and 



is the fermion number operator. Under- 



standing the properties of this model is a major challenge 
in the theoretical study of high temperature supercon- 
ductivity. In spite of the simple form of TL, solving the 
t-J model is non-trivial due to the strong interaction of 
the fermion objects. Mean-field theory solutions to the 
model often involve a fictitious statistical flux [Q. The 
flux pattern in these "flux phases" can be uniform or 
staggered. In particular, it has been proposed that the 
staggered- flux phase in the t-J model might lower the 
energy of the system and became the ground state at pa- 
rameters relevant to experimental systems ||. However, 
being an abstract mathematical quantity, it is difficult to 
find a suitable signature for the flux. Consequently inde- 
pendent confirmation of the existence of flux phases by 
numerical or experimental studies is difficult. A previous 
attempt to search for flux phases using exact diagonal- 
ization on small clusters has been either negative or 
inconclusive. 

As pointed out by Ivanov, Lee and Wen || , in the case 
of the doped t-J model a signature for the staggered-fiux 
phase can be found in the current-current correlation. 
Using a Gutzwiller-projected d-wave pairing wavefunc- 
tion, they found such a pattern in the current correla- 
tion, namely, the hole current goes around the elemen- 
tary square plaquettes in the counterclockwise and clock- 
wise directions, suggestive of positive and negative fluxes 
through the plaquettes in a staggered manner. They ex- 



plained this observation by showing that the Gutzwiller- 
projected d-wave pairing wavefunction is equivalent to 
the SU(2) projected staggered-flux wavefunction. Al- 
though this result shows that the concept of a staggered- 
flux phase is relevant in the projected d-wave pairing 
wavefunction, it does not answer the question whether 
this phase exists in the t-J model. The origin of the 
staggered-flux pattern in their study is in the staggered- 
flux phase of the SU (2) projected wavefunction. It will 
be very interesting to check whether such a pattern in 
the current correlation exists using an approach that is 
independent of the mean-field theory. 

Motivated by this result, we search for signatures for 
the staggered-flux phase in the t-J model using exact di- 
agonalization. This method is free from any analytical 
or numerical approximations, and thus can serve as an 
independent test for the result in Ref. ||. One serious 
drawback of exact diagonalization is that it is subjected 
to finite-size effect. This is particularly serious when we 
study correlation functions where the size and boundary 
conditions can have significant effects. To overcome the 
finite-size effect as much as possible, we use the largest 
cluster on which the t-J model is currently solvable by 
exact diagonalization, namely, the 32-site cluster with 
periodic boundary conditions. In order to study the cur- 
rent correlation, we need to dope the system with at least 
two holes. Such a two-hole t-J model on a 32-site clus- 
ter has recently been solved using exact diagonalization 
HQ . In this study we use the same lowest energy d-wave 
states at different J/t as in Ref. 0. They have d x 2_ y 2 
symmetry and are the respective ground states when J/t 
is not too small. Our primary interest is in the case 
J/t = 0.3 which is believed to be an appropriate value to 
describe the doped cuprates. With this parameter, the 
d-wave state is a weak two-hole bound state — in the 
sense that it has a negative but small two-hole binding 
energy E^/t = —0.05146. The root-mean-square separa- 



tion of the two holes is 



= 2.0587. Note that this 



value does not imply that the holes are tightly bound. 
Due to the finite system size, two uncorrelated holes on 
this cluster will have ^/(r 5 ) = 2.3827. 
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Once the desired wavefunction is found, it is straight- 
forward to evaluate the current correlation (jkijmn), 
where the current on a bond linking the sites k and I 
is defined by 



jkl = H^la^ ~ 4*Cko)- (2) 

Following the notation of Ref . || , we display the current 
correlation divided by the hole concentration, (jkijmn) /x. 
The result for J/t = 0.3 is shown in Fig. [|. It is to be 
compared with Fig. 1 of Ref. || . We note the very good 
agreement between these two results: not only all (ex- 
cept one) directions of the correlation agree, but also the 
corresponding numbers agree to the same order of mag- 
nitude. Given that the two results are obtained by very 
different approaches, their agreement is very surprising. 
The only exception to the agreement in the direction of 
the correlation is the one on the opposite side of the same 
square plaquette as the reference bond. This can be ex- 
plained by a feature in the hole-hole correlation function, 



C hh (r) = ((l-n r )(l-n )), 



(3) 



of the d- wave state that we use. While ChhO-) > Chh(V^) 
in the projected <i-wave pairing state the opposite is 
true in our d-wave state. In fact the holes in our rf-wave 
state have the largest probability to be at a distance \/2 
apart When the holes are at the opposite vertices 
of the same plaquette, the hole current along the two 
opposite sides of the plaquette should be along the same 
direction. 



vorticity V(r) of a square plaquette centered at r by sum- 
ming up the currents around it in the counterclockwise 
direction. The vorticity-vorticity correlation divided by 
the hole concentration, 



Cvv(r) = (V(r)V(0))/x, 



(4) 



is shown in Fig. ||(a). Note that Cyy can be computed 
from the current correlation, i.e., Fig. |^(a) is just another 
way of presenting the same data in Figs. |l|. Fig. ||(a) is 
to be compared to Fig. 1(b) of Ref. M. They all have 
the same staggered pattern. Note that this staggered 
pattern has a n phase shift, i.e., the sign of (V(r)V(O)) 
is ( — iy*+ r y+ 1 . This 7r phase shift was explained within 
the SU(2) picture in Ref. j|] by assuming the pairing of 
holes of opposite vorticities, i.e., circulating in opposite 
directions. 
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FIG. 2. The vorticity correlation Cvv(r) at (a) J/t = 0.3, 
and (b) J/t — 0.8. The reference plaquette at r = is indi- 
cated by a cross inside it. Only part of the cluster is shown. 
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FIG. 1. The current correlation (jkijmn) / 'x in the two-hole 
t-J model at J/t = 0.3. x = 1/16 is the hole concentration. 
The reference bond is indicated by an arrow with a circle. 
Arrows on other bonds point along the directions of positive 
correlation. Those arrows without numbers attached can be 
related by symmetry to those that have. 

In order to display the staggered pattern of the current 
correlation more clearly, we follow Ref. |j| and define the 



Since the hole current results from the motion of the 
holes, it is intuitively clear that the correlation between 
the current in different bonds must be related to the 
correlation between the holes in the corresponding lo- 
cations. Therefore, in addition to the above weak bound 
state at J/t = 0.3, we also study the wavefunction 
with the same symmetry but at J/t — 0.8 where the 
holes are more tightly bound with E^/t = —0.44423 and 
\J (r 2 ) = 1.5120. The purpose is to study how the cur- 
rent correlation changes when the binding of the holes is 
stronger. In Fig. ^(b) we show the vorticity correlation 
at J/t — 0.8. We note that it has the same staggered 
with a 7T phase shift pattern as that at J/t = 0.3. Since 
the holes are more tightly bound in the present case, we 
expect the current correlation, and therefore the vortic- 
ity correlation, to have a shorter range. This is indeed 
the case when we compare it with Fig. ||(a). Cyv(r) 
decreases much faster in the present case. 

In Ref. H it was observed that in the projected d- 
wave pairing state the hole-hole and vorticity-vorticity 
correlations decay with the same power law. In Fig. || we 
plot these two quantities in our d-wave states at J/t — 
0.3 and 0.8. As noted before Chh(r) depends on 

whether and r are on the same sublattice, and if so, 
whether r is an even (r x and r y are even) or odd (r x and 
r y are odd) site. Those on the same sublattice and odd 
sites are enhanced, while those on the same sublattice 



2 



and even sites are suppressed relative to those on different 
sublattices. If we assume that the hole correlation decays 
with a power law, Chh(r) ~ 1/?" Q , then the best fit gives 
a f» 1.15 < 2 at J/t — 0.3, making it an unbound state 
when the system size goes to infinity. Being a tightly 
bound state, the hole correlation at J/t = 0.8 decays 
much faster than that at J/t = 0.3 — the corresponding 
a « 3.04 is significantly larger than 2. The same trend 
is found for the vorticity correlation. Nevertheless, with 
the limited number of data points we are not able to say 
conclusively whether they decay with the same law, nor 
whether they decay with power laws at all. 
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FIG. 3. The hole correlation Chh{r) and the vorticity cor- 
relation |CVv(r)| of the d-wave states at J/t = 0.3 and 0.8. 
A (j4 s ) indicates that the pair of holes at and r are on the 
same sublattice and r is an odd (even) site. B indicates that 
they are on opposite sublattices. 

The above results show that a staggered-flux pattern 
exists in our exact cZ-wave state where the hole correla- 
tion decays with distance. In Ref. {m such a pattern is 
interpreted as the binding of holes of opposite vorticities 
in the SU (2) language. It is therefore interesting to see, 
using our exact numerical results, if such a pattern ex- 
ists in a state where the holes repel. For this purpose 
we use a low-lying excited state of the same two-hole 
system at J/t = 0.3. It is a doubly degenerate p-wave 
state with total momentum (-7T, tt) and is a spin singlet 
fiof . Fig. |i| shows the hole correlation of the p y -wave 
state. In contrary to the d-wave ground state, the holes 
in this state repel each other — Chh increases with r for 
small r. The mean-square separation between the holes is 
y/ (r 2 ) = 2.5376. Fig. || shows the current correlations in 
the p y -wave state. We observe that the magnitudes of the 
correlations are small even at small distances, and there 
is definitely no staggered-flux pattern. The same data 
plotted as vorticity correlation in Fig. |^ again confirms 
these observations. This is consistent with the interpre- 
tation that the staggered-flux pattern results from the 
pairing of the holes. 



FIG. 4. The hole correlation Chhij) of the p^-wave state 
at J/t = 0.3. Due to the lack of the rotational symmetry C4 
(i.e. a rotation by 7r/2) in the wavefunction, there may be two 
inequivalent points at the same r. Their values are indicated 
by filled circles. The solid line joins their mean values. 



FIG. 5. Same as Fig. |l| but for the p y -wave state. Due to 
the lack of the rotational symmetry C4, there is an inequiva- 
lent correlation where the reference bond is a vertical one. It 
is not shown for the sake of simplicity. 



FIG. 6. The vorticity correlation Cvvir) of the p y -wave 
state. 
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We now return to the d-wave state at J/t — 0.3 and 
search for another signature of the staggered-flux phase. 
It has been suggested that a physical interpretation of 
the statistical flux in the mean-field theory can be found 
in the spin chirality operator , 



X4 = 714 Si • (S 2 x S 3 ), 



(5) 



where 1, 2, 3, 4 are the vertices of a square plaque- 
tte in the counterclockwise direction. In Fig. ^ we plot 
the correlation C xx (r) = (xrXo) f° r the- d-wave state at 
J/t = 0.3. It is obvious that this correlation exhibits a 
short-range staggered pattern with a ir phase shift, but 
it decays very rapidly with distance. Since the chirality 
is related to a magnetic flux, which may in turn affect 
the movement of the holes, it is natural to associate the 
chiral order to the holes. To demonstrate the connec- 
tion between the chirality and the hole current, we follow 
Rcf. jl2| and define the spin chirality operator around a 
hole, 



x'l 



(l-n 4 )S 1 -(S 2 x S 3 ). 



(6) 



The chiral order set up by a hole is displayed in the cor- 
relation C x h x (r) = (XrXo) m Fig. 0. It has a staggered 
pattern with a tt phase shift. Compared to C xx (r), it 
decays slower, but has strong fluctuation as a function of 
r. Finally, we plot the correlation between the chirality 
around the two holes, C x h x h (r) = (XrXo)- decay with 
distance is again much slower than in the case of C xx (r) . 
These correlations show that the chiral order must be tied 
to the holes. Unfortunately, with the limited data we are 
not able to deduce conclusively whether long-range chi- 
ral order exists in the current system. We also note that 
similar results on the chirality correlations were obtained 
with the projected e?-wave pairing wavefunction jlifl . 
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FIG. 7. Correlations of the chirality operators of the d-wave 
state at J/t — 0.3 in the same logarithmic scale. Empty 
and solid circles represent positive and negative values re- 
spectively. 
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current and chirality correlations in the lowest-energy d- 
wave state of the two-hole t-J model on a 32-site cluster 
with periodic boundary conditions. The current corre- 
lation clearly shows a staggered-flux pattern. Note that 
in the d-wave states at different J/t, the holes are mu- 
tually attractive — Chh(f) decays with distance. In 
this case the staggered correlation in the vorticity de- 
cays faster when the attraction between the holes be- 
comes stronger. In the p-wave state where the holes re- 
pel {Chh{r) increases with distance at small distances), 
no such staggered-flux pattern is found. Chirality corre- 
lations also show a similar staggered pattern. Although 
the numbers are small, they show that any long-range 
pattern in the chirality correlation must be related to 
the holes. Our results on the current correlation are con- 
sistent with a recent SU (2) study by Ivanov, Lee and 
Wen ||. Being a completely independent approach, the 
fact that our unbiased numerical result agrees with theirs 
gives strong support to the notion that the staggered-flux 
phase is a key ingredient of the ground state in the t-J 
model, at least in the low doping regime. 
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In conclusion, we have looked for signatures for the 
staggered-flux phase in the t-J model by calculating the 



4 



